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ABSTRACT 

Action functionals describing relativistic perfect fluids are presented. Two of these 
actions apply to fluids whose equations of state are specifled by giving the fluid 
energy density as a function of particle number density and entropy per particle. 
Other actions apply to fluids whose equations of state are specifled in terms of 
other choices of dependent and independent fluid variables. Particular cases include 
actions for isentropic fluids and pressureless dust. The canonical Hamiltonian forms 
of these actions are derived, symmetries and conserved charges are identifled, and 
the boundary value and initial value problems are discussed. As in previous works 
on perfect fluid actions, the action functionals considered here depend on certain 
Lagrange multipliers and Lagrangian coordinate flelds. Particular attention is paid 
to the interpretations of these variables and to their relationships to the physical 
properties of the fluid. 



1. Introduction 



This paper contains a description of several action functionals whose extrema 
are the classical histories of a relativistic perfect fluid. The possible benefits of 
these actions are difficult to judge, but experience has shown generally that action 
principles are a powerful and conceptually elegant means of specifying a dynamical 
system [1]. They concisely encode the classical equations of motion, provide a 
direct and simple way of relating symmetries and conserved charges, and serve as 
the starting point for a canonical Hamiltonian analysis and for Hamilton-Jacobi 
theory. For a quantum mechanical description of a system, the existence of an 
action functional appears to be a necessity [2]. 

Perfect fluids are described locally by various thermodynamical variables. In 
the notation of Misner, Thorne, and Wheeler [3], they are 



n = particle number density , (l.lo) 

p — energy density , (1-1^) 
p = pressure , 

T = temperature , (l.ld) 

s = entropy per particle . (l-l^) 



These variables are spacetime scalar fields whose values represent measurements 
made in the rest frame of the fluid. The fluid motion can be characterized by its 
unit four- velocity vector field U"'.* 

Previous works on relativistic and nonrelativistic perfect fiuid actions [4-13] 
have revealed a number of common features that any such action principle must 
possess. In particular, it is known that no perfect fiuid action can be constructed 
solely from the variables (1.1) and U°' unless the variations among those variables are 
constrained [5]. Two of the required constraints are particle number conservation 
{nU"')-a = and the absence of entropy exchange between neighboring fiow lines 
{nsU°')-a = 0. The remaining constraint is that the fiuid fiow lines should be fixed 
on the boundaries of spacetime. One way of enforcing this constraint is to vary 

* Spacetime tensor indices are denoted by a, 6, etc., the spacetime metric is 
Qab, and spacetime covariant derivatives are denoted by a semicolon. The sign 
conventions of reference [3] are used thoughout. 
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the flow lines by Lie dragging along some vector field where vanishes on the 
boundaries. The flow line variation then induces variations in the four-velocity U"' 
with the result SU"- = (5^ + UbU"-)£^U^, where £^ is the Lie derivative along 

There is another method for handling the constraint that the flow lines should 
be fixed on the spacetime boundaries, suggested by the work of Lin [9, 10] . Instead of 
the four- velocity t/", the history of the fluid is characterized by a set of spacetime 
scalar flelds a^, A = 1, 2, 3, that are interpreted as Lagrangian coordinates for 
the fluid. That is, a^{x) serve as labels for the fluid, specifying which flow line 
passes through a given spacetime point x. A set of Lagrangian coordinates can be 
generated by choosing an arbitrary spacelike hypersurface and a coordinate system 
on that surface. Then each flow line is labeled by the coordinate value of the 
point where it intersects the hypersurface. By building an action functional using 
the Lagrangian coordinates a^, the fluid flow lines are held fixed on the spacetime 
boundaries by simply fixing on the boundaries. 

The particle number and entropy exchange constraints can be incorporated 
directly into the action via Lagrange multipliers. An action S is presented in section 
2 that makes use of such Lagrange multipliers along with the Lagrangian coordinates 
to enforce the constraints. This action describes a perfect fluid whose equation 
of state is specified by giving the energy density p as a function of number density 
n and entropy per particle s. The action S incorporates various features found in 
the action functionals discussed in references [5, 8, 14-16], but to my knowledge it 
has not been written previously in precisely the form discussed here. In sections 
2 and 3 the Lagrangian and Hamiltonian forms of the action S are developed, the 
Lagrangian and Hamiltonian equations of motion are displayed explicitly, global 
symmetries and conserved charges are analyzed, and the initial and boundary value 
problems are discussed. 

The use of Lagrangian coordinates ct"* and Lagrange multipliers is sometimes 
criticized on the grounds that these extra variables are "unphysical" . One of the 
results of the present work is to show that the values of each of these variables are 
determined to within a global symmetry transformation by the physical properties 
of the fluid. The situation here is formally analogous to that of a free nonrelativistic 
point particle moving in flat space with cartesian coordinates. In that case, New- 
ton's law is = 0, where v is the particle's velocity and the dot denotes a time 
derivative. An action functional that yields this equation of motion is J dt{v ■ v/2), 
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where the variations in v are constrained to those having the form Sv = ^. Here, ^ is 
an infinitesimal vector that vanishes at the initial and final times. An unconstrained 
action is obtained by introducing a "Lagrangian" coordinate x and constructing the 
velocity from v = x. The new variable x is just the particle's coordinate location. 
The value of x at any given time specifies the physical location of the particle to 
within a symmetry translation or rotation of the cartesian coordinates. Similarly, 
the values of the Lagrangian coordinates specify the physical location of the 
fluid flow lines to within a symmetry transformation that amounts to a change of 
Lagrangian coordinates. The values of the Lagrange multipliers that appear in the 
fluid action S are also determined to within a global symmetry transformation by 
certain physical properties of the fluid. (The action S is the analogue of the free 
particle action J dt{v ■ x — v ■ v/2) in the sense that it is a functional of both the 
four-velocity U°' and the Lagrangian coordinates a^.) 

One of the fluid equations of motion obtained from the action S relates the 
fluid four-velocity Ua to the Lagrange multipliers, Lagrangian coordinates a^, and 
their gradients. Expressions of this type are often called velocity-potential repre- 
sentations or Clebsch [29] representations of the four-velocity. There have been 
numerous discussions in the literature concerning both the number of scalar flelds 
that are mathematically required for such a representation, and the number of scalar 
flelds that are physically natural for such a representation [4, 5, 10, 12, 15]. Section 
4 contains a further discussion of these issues. 

The Hamiltonian form of the perfect fluid action is a functional of the fluid 
number density, the fluid entropy density, the Lagrangian coordinates a'^, and their 
canonical conjugates. In this case the number and entropy densities are Euler- 
ian, that is, their values correspond to measurements made by observers at rest in 
space. The Eulerian densities arc related to the Lagrangian (or comoving) densities 
n and ns by a kinematical "gamma" or boost factor that in turn is determined by 
the local spatial velocity of the fluid. The fluid contributions to the Hamiltonian 
and momentum constraints are just the appropriate projections of the perfect fluid 
stress-energy-momentum tensor [17], and involve the spatial components Ui of the 
fluid four-velocity. These components Ui are explicitly expressed in terms of the 
canonical fluid variables and the Lagrangian particle number density n. The La- 
grangian number density n is itself a function of the canonical variables as implicitly 
determined by the equation relating the Lagrangian and Eulerian number densities. 
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A Hamiltonian formulation of perfect fluids that makes use of only "physical" 
fluid variables has been developed (see [18] and references therein) based on the 
Lie-Poisson brackets. However, Lie-Poisson brackets are not canonical, so there is 
no underlying symplectic structure and the usual "J {pq — H)" form of the action 
does not apply. On the other hand, the Lie-Poisson structure can be derived by a 
so-called Lagrangian to Eulerian map starting from a certain canonical Hamiltonian 
structure that involves the Lagrangian coordinates [18]. This canonical Hamilto- 
nian formulation of perfect fluids is derived in section 5, starting from an associated 
action S that is a functional of the Lagrangian coordinates only. In this action, 
the particle number and entropy exchange constraints are enforced by restricting 
nU^ and s to depend on spacetime only through certain combinations of the Lagan- 
gian coordinates. In particular, s is given by a function s{a), and nil"' is given by 
—'ni23e°'^'^'^o^,bO?,cOi^,di where 77123 is a function of a^. {e"'^'^'^ is obtained by raising 
indices on the totally antisymmetric spacetime volume form eabcd-) The action S, 
like S, describes a perfect fluid whose equation of state is specifled by giving the 
energy density p as a function of number density n and entropy per particle s. 

In addition to the local thermodynamical variables (1.1), it is also convenient 
to deflne [3] 

fi := ^ ^ = chemical potential , (1.2a) 
n 

a := — —Ts = physical free energy , (1.2b) 

n 

f := Ts = chemical free energy . (1.2c) 

n 

The chemical potential // is the energy per particle required to inject a small amount 
of fluid into a fluid sample, keeping the sample volume and the entropy per parti- 
cle s constant [3]. Similarly, the physical free energy a is the injection energy at 
constant number density and constant total entropy, while the chemical free energy 
/ is the injection energy at constant volume and constant total entropy. The ther- 
modynamical variables (1.1-2) are related by the local expression of the first law of 
thermodynamics, namely [3] 

dp = jj, dn + nT ds . (1-3) 

This relationship shows that an equation of state for a perfect fluid can be specified 
by giving the function p{n,s), the energy density as a function of number density 
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and entropy per particle. The first law also can be written as 

dp = n dfj, — nT ds (1-4) 

which naturally suggests an equation of state of the formp(//, s). Another possibility 
is 

d{na) — f dn — ns dT (1-5) 

which suggests an equation of state a(n, T). In section 6, various action functionals 
are presented for perfect fluids specified by the equations of state s) and a(n, T). 
As special cases, actions for isentropic fluids and pressureless dust are obtained. 

The perfect fluid action functionals developed by Taub [6, 7] and Schutz [15, 
19, 20] have been applied to the analysis of the stability of stellar models [7, 21-23], 
and it is hoped that the actions discussed here can be put to similar use. Schutz's 
action also has been used to quantize the combined gravity-fluid system [24, 25], 
where one of the Lagrange multipliers is chosen as a time coordinate. Along these 
same lines, the action for a pressureless dust, discussed in section 6, has been used in 
an investigation of the problem of time in canonical quantum gravity [26] . Another 
application of perfect fluid actions is found in the work of Gibbons and Hawking 
[27] (also see [28] and references therein). There, the action is used to obtain a 
thermodynamical potential that characterizes the global thermal properties of a 
gravitating fluid system. 

2. Action S with equation of state p(n, s) 
2.1 Action and equations of motion 

By definition, the stress-energy-momentum tensor for a perfect fluid has the 

form 

rj^ab ^ pjjajjb ^ ^(^a6 ^ ^a^6) ^ (2.1) 

where U"" is the unit four-velocity of the fluid. The equations of motion for a 
perfect fluid consist of the stress tensor equation of motion, namely T"^;;, = 0, and 
the equation {nil"') -a = expressing conservation of particle number [3]. The action 
functional presented in this section incorporates the stress tensor (2.1), the required 
equations of motion, and the first law of thermodynamics (1.3) for a perfect fluid 
with equation of state p(n, s). 
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The perfect fluid action 5 is a functional of a spacetime contravariant vector 
density J"' that is interpreted as the densitized particle number flux vector -y/^nt/". 
That is, the fluid four- velocity is deflned by 

U" := J'*/! J| , (2.2) 

where 



\J\ := V^J^^g^ (2.3) 

is the magnitude of J", and the particle number density is given by 

n := \J\/V^ . (2.4) 

This action is also a functional of the spacetime metric gab, the entropy per particle 
s, the Lagrangian coordinates a^, and spacetime scalars denoted by ip, 6, and Pa- 
(The indices A, B take the values 1, 2, 3.) In terms of an arbitrary equation of 
state p{n,s), the action reads 

+ J''i<f,a + se,a + P^a\a)} ■ (2.5) 

The stress-energy-momentum tensor derived from this action is 



.aH 2 SS 



and has the perfect fluid form (2.1) with energy density p and pressure defined by 

p:=„|^-p. (2.7) 

This expression for pressure agrees with the relationship implied by the first law of 
thermodynamics (1.3). 

The fluid equations of motion derived from the action (2.5) are as follows: 

0=^ =//l7« + + , (2.8a) 

0=g = -iV, (2.86) 
=^ = -(sJ"),. , (2.8c) 

0=^ = -(;3./°),a, (2.8e) 
0=^ = a^.^°. (2.8/) 



Solutions to these equations extremize the action S with respect to variations that 
leave (f, $, and a"^ fixed on the spacetime boundaries. According to equations 
(2.8b,c), the fields (p and 9 serve as Lagrange multipliers for the particle number 
conservation constraint J^^a = and the entropy exchange constraint {sJ°')^a = 
0, respectively. Equation (2.8f) shows that /?a are Lagrange multipliers for the 
constraints a^^aJ" = that restrict the fluid four-velocity vector to be directed 
along the flow lines = constant. 

As discussed in the introduction, the fields a^{x) are interpreted as Lagrangian 
coordinates for the fluid and serve as labels that specify which flow line passes 
through a given spacetime point x. A set of Lagrangian coordinates can be generated 
by choosing an arbitrary spacelike hypersurface and specifying a coordinate system 
CK^ on that surface. Then each flow line is labeled by the coordinate value of 
the point where it intersects the hypersurface. It is useful to view this arbitrary 
spacelike hypersurface as an abstract "fluid space" whose points represent the fluid 
flow lines, and to view the Lagrangian coordinates as a coordinate system on 
the fluid space. (The concept of "fluid space" is also discussed in reference [14], 
where it is called the "matter space".) So conceptually, the fluid space is the space 
of Lagrangian coordinate labels, and as a manifold it is isomorphic to any spacelike 
hypersurface. It may be impossible to cover the fluid space with a single coordinate 
chart, so the Lagrangian coordinates generally must be deflned in open subsets 
of the fluid space. With the above interpretation, it is assumed that the flelds 
constitute a good set of Lagrangian coordinates, in the sense that within the 
appropriate open subsets each flow line carries a unique label and the gradient 
is nonvanishing. 

The significance of the scalar field 9 is revealed by comparing the equation 
of motion (2.8d) with the first law of thermodynamics (1.3). This leads to the 
identification 

T = 9,aU^ = -^ (2.9) 

n OS 

for the fluid temperature T, and shows that 9 is the "thermasy" discussed by van 
Dantzig [30]. Kijowski et al. [14] have proposed a physical interpretation for ther- 
masy, relating it to the difi^erence between proper time along the fluid flow lines and 
proper time along the history of a typical fluid particle that executes chaotic motion 
about the flow lines. The field ip plays a role that is mathematically analogous to 
9, as seen by contracting equation (2.8a) with t/" and using equations (2.8f) and 
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(2.9) to obtain 



(2.10) 



Thus, is a "potential" for the temperature T and (pis a "potential" for the chemical 
free energy /. The interpretations of the fields 9, (p, and Pa sixe discussed further 
in sections 3 and 4. 

As stated above, the perfect fluid equations of motion (2.8) should imply 
particle number conservation and the stress tensor equation of motion T"'^-b ~ 0. 
The first of these, particle number conservation, is expressed explicitly in equation 
(2.8b). For the stress tensor equation of motion, first consider its projection along 
the fluid flow lines. This yields [3] 

= C/,T«% = -|^s,„C/" , (2.11) 

where particle number conservation has been used. This equation combines parti- 
cle number conservation (2.8b) and the entropy exchange constraint (2.8c), and is 
indeed implied by the equations of motion. This result shows that the fluid flow 
is locally adiabatic, that is, the entropy per particle along the fluid flow lines is 
conserved. 

The projection of the stress tensor equation of motion orthogonal to the flow 
lines gives the Euler equation, relating the fluid acceleration to the gradient of 
pressure [3]: 

0^{9ab + UaUb)T'-., 

= (P + p)Ua;bU'' + {Si + UaU^)p,b . (2.12) 

Using expression (2.7) for p, the Euler equation can be written as 

O^I^Ua.,bU'+{Si + UaU')-p,b 

= 2{fxU,a).MU' - {Si + UaU')^^s,t , (2.13) 
9 



where the square brackets denote antisymmetrization. Now combine equation (2.13) 
with the entropy conservation equation (2.11) to obtain 

2V^a-MU'' = Ts,a , (2.14) 

where 

Va := fiUa (2.15) 

and the definition (2.9) for temperature T has been used. The vector Va, sometimes 
called the Taub current, plays an important role in relativistic fluid dynamics, es- 
pecially in the description of circulation and vorticity [31]. From the interpretation 
of fi as an injection energy, V"" is identified with the four-momentum per particle of 
a small amount of fluid to be injected in a larger sample of fluid without changing 
the total fluid volume or the entropy per particle. 

It remains to be verified that the equations of motion (2.8) imply the Euler 
equation (2.12). For this purpose, it suffices to show that the equations of motion 
imply equation (2.14), since that equation and the Euler equation are related by 
use of the equations of motion. Accordingly, compute 2V[a;b]U^ using the expression 
for Va from equation (2.8a) to obtain 

^V^a^U" = -2(^,[a + se,^a + P^a^ ,la)MU'' • (2-16) 

The equations of motion imply that the terms involving gradients of s, and 
along the flow lines all vanish. Only the term s^ad,bU^ remains on the right-hand- 
side of equation (2.16). According to equation (2.9), this term equals Ts^a so that 
the equations of motion indeed imply equation (2.14). The conclusion is that the 
equations of motion (2.8) derived from the action (2.5) imply the relevant perfect 
fluid equations of motion, including the Euler equation. 

The "on shell" action, that is, the value of the action (2.5) when the equations 
of motion (2.8) hold, is just the proper volume integral of the pressure: 



^(on shefl) = 




(2.17) 



Of course, the addition of surface integrals to the action will change its on-shell 
value without affecting the equations of motion. In particular, the surface integral 
a can be added to S, which amounts to replacing the term ip^aJ"' 
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by the term —ipj^-^a- Likewise, adding the surface integral — J d'^x{9sJ°')^a to 5' 
results in replacing sd^aJ"' by —9{sJ°')^a- These surface integrals have on-shell 
values — J d^x^^/—gnf and — J d^x^J—gnTs, respectively. Thus, by adding these 
surface integrals in various combinations, action functionals can be obtained whose 
on-shell values are, for example, the proper volume integrals of —p or —na. 

There are a number of alternative perfect fluid actions that differ from the 
action (2.5) in the placement of derivatives in each term. In some cases, such as 
changing ^p^aJ"' to —ipJ°'^ai this amounts to the simple addition to the action of a 
boundary term as discussed above. A valid action functional is also obtained by re- 
placing sd^aJ"' with —ds^aJ"'- This change is brought about by defining (p := (p' — s6 
then dropping the prime on the new field {(p' — > <^). In place of equation (2.10), 
(p then satisfies /j, = (p^aU"' and is a "potential" for the chemical potential. Nev- 
ertheless, the resulting action yields the required equations of motion and stress 
tensor. The action (2.5) has an advantage over its alternatives in that the momenta 
conjugate to (p and $ that naturally follow from that action have clear physical in- 
terpretations as the particle number density and entropy density seen by (Eulerian) 
observers at rest in space. 

2.2 Symmetries of the action S 

The perfect fluid action (2.5) is invariant under the inflnitesimal transforma- 
tions* 

5(p = eF- es{dF/ds) - tl3j,{dF/dl3j,) (2.18a) 

59 = e{dF/ds) (2.186) 

5s = (2.18c) 

5a^ = e{dF/dp^) (2.18d) 

5p^ = -e{dF/da^) , (2.18e) 

along with 5 J" = and 5 gab = 0, where e is an inflnitesimal parameter and F is 
an arbitrary function of the Lagrangian coordinates a^, the entropy per particle 
s, and These invariances, one for each function F{a,P,s), constitute global 
symmetries of the theory, symmetries in which the fleld transformations are flxed 



* These symmetries are derived in the Appendix. 
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for all time. Such a symmetry is to be distinguished from a gauge symmetry, in 
which the transformation can vary from one instant of time to the next. Observe 
that the particular combination (fi,a + s$^a + PAOi^,a that appears in the action (2.5), 
and equals —jjXJa when the equation of motion (2.8a) holds, is invariant under 
the transformations (2.18). Also note that the function F{a,P,s) has vanishing 
gradient F^glJ"' = along the fluid flow lines, according to the equations of motion 
(2.8). Thus, the symmetries (2.18) produce variations in ^, a^, and (3 a that 
are constant along the flow lines. Such variations do not affect the temperature 
T = 0,aU°'i chemical free energy / = (p^aU"', or the constancy of and along 
the flow lines. 

The physical meaning of the symmetries (2.18) can be explored by consider- 
ing special choices for the function F{a,P,s). If F is a function of only, the 
transformations (2.18) reduce to 

Sip = eF (2.19a) 
5j3A = -e(dF/da^) , (2.196) 

with all other fields unchanged. Since the Lagrangian coordinates uniquely label 
the flow lines, at least within open subsets of the fluid space, the change (2.19a) in 
ip along one flow line is independent of the changes in (p along the other flow lines 
(subject to the requirement of continuity in (p). Therefore, the symmetries (2.19) 
can be understood as arising from the freedom to shift the value of (p along each 
flow line by a constant amount. In particular, any solution to the fluid equations 
of motion can be transformed via equations (2.19) into a solution with, say, (p = 
on any given spacelike hypersurface. 

If the function F has the form s F{a) for some function F of a^, the trans- 
formations (2.18) reduce to 

Se = eF (2.20a) 
5/3^ = -es{dF/da^) , (2.206) 

with the other fields unchanged. These symmetries can be used to shift the value 
of the thermasy 6 along each flow line by a constant amount. Thus, any solution 
of the equations of motion can be transformed into a solution with ^ = on any 
given spacelike hypersurface. 
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Now consider functions F of the form P^F^{a) where is a set of functions 
of a^. In this case, the transformations (2.18) reduce to 



eF^ 



(2.21a) 
(2.216) 



5(3 J, = -e(3s{dF''/da'') . 



These symmetry transformations describe changes of coordinates a"^ in the fluid 
space, where F^ is viewed as a vector in the fluid space. Moreover, observe that the 
equations of motion (2.8b,e) imply the constancy of 13a along the fluid flow lines, 
so that 13a can be expressed as a function of the Lagrangian coordinates a^. Then 
according to its transformation (2.21b) under changes of coordinates (2.21a), 13a 
can be viewed as the components of the covariant vector (or one-form) j3 :— (3Ada^ 
in the fluid space. 

The conserved Noether currents associated with the symmetries (2.18) are 
obtained from the general variation of the action (2.5), which is 



The action S is invariant under the symmetry transformations (2.18), so when the 
equations of motion hold, = 5S becomes 



Therefore FJ" are the conserved currents that satisfy {FJ"')^a = by virtue of the 
equations of motion (2.8). If space is closed, or if suitable boundary conditions are 
imposed at spatial infinity, the volume integral (2.23) can be written as a surface 
integral Q[F] evaluated on a final spacelike hypersurface, minus the same surface 
integral evaluated on an initial spacelike hypersurface. That surface integral is the 
Noether charge 



and equation (2.23) just expresses the conservation of Qf-F], that is, the indepen- 
dence of Q[F] on the choice of hypersurface E. In equation (2.24), h is the de- 
terminant of the three-metric on E, and Ua is the future pointing unit normal to 



5S = (terms giving the equations of motion) 




(2.22) 




(2.23) 




(2.24) 
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Insight into the meaning of the conserved charges Q[F] is obtained by consid- 
ering again a special choice for F{a,P,s). Let F be a function of only, and in 
particular choose F equal to unity for all in some ball B in the fluid space, and 
equal to zero outside B. Then the charge (2.24) becomes a proper volume integral 
over the subspace of the hypersurface S that contains flow lines in B, with integrand 
n{—naU°'). The factor —UaU"' is the relativistic "gamma factor" characterizing a 
boost from the Lagrangian observers with four-velocity U"' to the Eulerian observers 
with four-velocity [32]. Thus, n{—naU°') is the particle number density as seen 
by the Eulerian observers. Then conservation of the charge (2.24) expresses the 
conservation of particle number within a flow tube defined by the bundle of flow 
lines contained in the ball B. 

For general functions F{a,P,s), the conserved charges (2.24) can be given a 
physical interpretation in terms of the case discussed above. Since and s are 
constant along the flow lines, and because are unique flow line labels, the fields 
Pa and s can be expressed as functions of the spacetime point x through the com- 
bination a^{x). Thus, a general function F can be viewed as having functional 
dependence F{a, P{a), s{a)). Then the charge Q[F] is the number of particles con- 
tained in a distribution of flow lines, with each flow line weighted by the distribution 
function F{a,P{a),s{a)). For example, choose F equal to sF, where F equals unity 
for all in a ball 5, and zero outside B. In this case the conserved charge Q[F] 
equals the number of particles contained in the flow lines included in S, weighted 
by the entropy per particle s. Conservation of Q[F] expresses the conservation of 
total entropy within the flow tube deflned by the bundle of flow lines contained in 
B. 

The conserved charges (2.24) must be distinguished from quantities that are 
conserved in a more restricted sense, such as the fluid circulation C. Circulation is 
deflned as the integral of the Taub current one-form Va = jJiUa around a spacelike 
loop [31]. C is conserved in the sense that its value is a constant for particular 
families of loops, namely those obtained by evolving an initial loop along the fluid 
flow lines by an amount proportional to the thermasy 9 [33] . (If the loop is evolved 
along the flow lines by an amount proportional to proper time, changes in C are 
determined by the gradient of the entropy per particle [31].) On the other hand, a 
conserved charge Q[F] has the same value for any spacelike hypersurface, not just 
for a particular family of hypersurfaces. See reference [34] for a discussion of the 
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various types of conservation laws that arise in fluid mechanics. 



3. Hamiltonian form of the action S 

3.1 Action and Hamiltonian 

The perfect fluid action (2.5) can be written in Hamiltonian form by intro- 
ducing first a space-time decomposition of the fields. Accordingly, let t denote a 
scalar function on spacetime, whose gradient is nonzero and (future pointing) time- 
like. The t = constant surfaces foliate spacetime into spacelike hypersurfaces with 
unit normal = —Nt^ai where N := (— t,afi'"''^,fe)~^^^ defines the lapse function. 
Also introduce a time vector field t"" such that t^at"' = 1 and the projection tensor 
^& + n"nb onto the leaves of the foliation. With the shift vector defined by 
:= h^t^, the time vector field can be written as = Nn"" + N"-. 
Using the above relationships, the vector density J" is decomposed as 

= {hjy + {t^- - iv")n , (3.1) 

where the definitions (hJ)"" := h^J^ and 11 := t^aJ"" are used. Because J" = 
^/—gnU'^ is the densitized particle number fiux vector, this latter definition becomes 

n = y/hni-UaU") , (3.2) 

where ^y—g = N\/h has been used. Recall that —UaU^ is the relativistic "gamma 
factor" characterizing a boost from the Lagrangian observers with four-velocity 
to the Eulerian observers with four-velocity n". Then 11 is recognized as the 
spatially densitized Eulerian particle number density; that is, H/^/h is the number 
density of fluid particles as seen by the Eulerian observers who are at rest in the 
t — constant hypersurfaces. 

With the decomposition (3.1), the action (2.5) becomes 

s^J d^x[-^p{\j\/^, s) + m«(^,„ + se,a + pAa\a) 

+ {{hJr-UN'^){<fi,a + se,a + P^a^,a)} ■ (3.3) 

15 



Now tie the spacetime coordinates to the fohation by choosing t as the "timehke" 
coordinate, and choosing spatial coordinates such that {d/dtY = t°'. In these 
coordinates, the spacetime metric has the ADM form [35] 

gabdx^dx^ = -N'^dt^ + % [dx' + N'dt){dx^ + N^dt) , (3.4) 

where a;* are the spatial coordinates on the t = constant hypersurfaces, and hij are 
the spatial components of the spacetime tensor hat- Note that derivatives along 
are just ordinary t derivatives, and will be denoted by a dot. Since the contravariant 
tensors N"- and {hJ)°' have vanishing t-components, the fluid action (3.3) can be 
written as 

S=J dtd^x^n{(p + se + P^a^) - N^/hp{\J\/N^/h,s) 

+ {{hjy-UN'){if,i + se,i + p^a^,i} , (3.5) 

where 



\J\ = 



{Nnf -{hjy{hj)iY^'' , (3.6) 



is the magnitude of J". (Spatial indices are raised and lowered by the spatial metric 
hij and its inverse h'^^ .) 

Since the fields with a time derivative appear linearly in the action (3.5), it 
is clear that 11 is canonically conjugate to (p, Us is canonically conjugate to 9, 
and UPa is canonically conjugate to a^. Only the fields {hjy have no canonical 
counterpart. The action can be cast into canonical form with respect to the variables 
{hjy by following the standard analysis due to Dirac [36]. The momenta conjugate 
to (hjy are constrained to vanish, and the preservation in time of these primary 
constraints is equivalent to the equation of motion produced by varying the action 
with respect to (/iJ)*; this is just the projection of the equation of motion (2.8a) 
onto the spacelike hypersurfaces, namely 

Ui = -{ip^i + s9,i + /3Aa^,i)/fi. (3.7) 

The complete set of constraints is second class, and their elimination is equivalent 
to replacing {hjy in the action (3.5) by the solution to the equation of motion (3.7). 

Equation (3.7) can be solved implicitly for (hjy as a function of 11, 6, s, 
a^, and hij in the following way. Observe that {hjy appears on the left-hand- 
side of that equation through the combination Ui = {hJ)i/\J\, with |J| given by 
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equation (3.6). The equation Ui = {hJ)i/\J\ can be solved for {hjy as a function 
of t/j, with the result 

NU 

(hjy = W . (3.8) 

Now consider Ui as shorthand notation for the right-hand-side of equation (3.7). 
This expression for Ui depends on the chemical potential /j, = dp/dn, which is a 
function of n and s as determined by the equation of state p{n,s). In turn, n is 
determined implicitly as a function of H, (p, $, s, P^, a^, and hij by the equation 

n = Vhn{l + Uih'WjY''^ . (3.9) 

This expression (3.9) relates the Lagrangian number density n to the Eulerian num- 
ber density 11/ \/h, and is obtained from equation (3.2) with the gamma factor writ- 
ten as -iiaU'' = (1 + Uih^Wjy^. The relationship -UaU" = (1 + Uih^^Ujf^ is 
proved by inserting the space-time split of the inverse metric g""^ into the identity 
— 1 = Uag"'^Ub and solving for (1 + UiK^Wj). To summarize, Ui as it appears in 
equations (3.8) and (3.9) should be viewed as shorthand notation for the right- 
hand-side of equation (3.7). Then equations (3.8) and (3.9) implicitly determine 
[hjy and the Lagrangian number density n as functions of the Eulerian number 
density H, the fields 6*, s a^, and the spatial metric hij. 

The Hamiltonian form of the action is now obtained by substituting expression 
(3.8) for {hjy into the action (3.5) and using equation (3.7). This yields 

S = j dtd^x^ip + {Iis)e + (n/3^)d^ - W""^ - AT^Tif , (3.10) 

where the fluid contributions to the Hamiltonian and momentum constraints are 

= Vh[p{l + U'Ui) + pU'Ui] , (3.11a) 
^ _^nt/i . (3.116) 

The corresponding perfect fluid Hamiltonian is 

^fiuid _ j d^x{Nn''""' + N'nf''"') . (3.12) 

In equations (3.10-12), n and Ui are considered to be functions of the variables H, 
(f, {Us), 9, (n/3^), a"^, and hij as determined by equations (3.7) and (3.9). 
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The fundamental Poisson brackets among the fluid variables are 

{ip{x),Il{x')} = d{x,x') , (3.13a) 
{e{x), {Ils){x')} = S{x,x') , (3.136) 
{a^{x), {Up^){x')} = SiS{x,x') . (3.13c) 

Observe that the fluid contributions (3.11) to the Hamiltonian and momentum con- 
straints are given by the projections TY""''^ = VhuaT'^'^ni, and TYf"''* = \/hnaT'^i 
of the stress tensor (2.1), as expected of a matter action with nonderivative cou- 
pling to gravity [17]. Correspondingly, the gravitational fleld contributions to the 
constraints are the same as for vacuum general relativity. The fluid contribution 
(3.11a) to the Hamiltonian constraint can be rewritten in various useful forms by 
using expression (3.9) to relate the Eulerian and Lagrangian number densities: 

= iM^/iVhn) - Vhp . (3.146) 



Also observe that the momentum constraint (3.11b) equals 



(3.15) 



which is the form dictated by the role of TY""''* as the generator of spatial diffeomor- 
phisms for the scalar fields ^p, 9, and and their conjugates. 

3.2 Canonical equations of motion and symmetries 

In order to compute the canonical fiuid equations of motion, first vary expres- 
sion (3.14b) for 7^*1"'"^ with respect to n, s, H, and hij. This yields 



-nil + U -U) +n{U -U) 



dn + Vh 



+ ^^mii + ]-\fh[-nii{l + U ■ U)-p]W^hij , 
\'hn 2 



Ss 
(3.16) 



where U ■ U := W^Ui and :— h^^Uj, and expression (3.9) is used. By varying 
equation (3.9) for H, Sn is given in terms of variations in the canonical variables by 



^Vh 



-l^{l + U -U) +n{U -U) 



+ Vhn(U ■ U)^—^Ss + ii(l + U- Uy/^SU - VhnWSUUi) 
onos 



6n - -Vhnp[{l + U ■ U)h'^ - WW]6hi 

(3.17) 
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Combining these results gives 

+ (1 + U- U)-^/^ [fSU + T5(ns) - ^] , (3.18) 

where Ui and n are functions of the canonical variables as determined by equations 
(3.7) and (3.9). In addition, the pressure p, temperature T, and chemical free 
energy / are functions of the canonical variables defined through the relations p = 
n{dp/dn) — p,T = [dp/ ds) /n, and / = (dp/dn) — Ts. 

The functional derivatives of the fluid Hamiltonian (3.12) follow from the 
variation (3.18) and the variation of Ti""''^ from equation (3.15). In canonical form 
the perfect fluid equations of motion are 

^^fluid 



N{f - ^7V,^)(1 + U ■ t/)-V2 + iVV,i , (3.19a) 



e = j^—y = N{T - we,i){i + u ■ u)-^/^ + N'e,i , (3.196) 

-1 

N{-Wa'',i){l + U ■ U)-^/^ + N'a^^i , (3.19c) 



^^fluid 



X tr fluid 

n=-^ ^-{NVhnU')i + {Nm)i , (3.19d) 

X tJ- fluid 

(UsY = = -{N^/hnsW),, + (iV^Hs),, , (3.19e) 

(Up^y = = -(NVhnp^W),, + (Nmp^),, . (3.19/) 

These are precisely equations (2.10), (2.9), (2.8f), (2.8b), (2.8c), and (2.8e), respec- 
tively, written in terms of the 3+1 decomposition (3.1), (3.8) for J". Also observe 
that the variation (3.18) yields 



J^flui, 



6hij 2 



--N^/h[ph'^ + {p + p)WW] , (3.20) 



which is the fluid contribution to the canonical equation of motion for the spatial 



metric hij. 



The global symmetries (2.18) described in section 2.2 appear in the canoni- 
cal formalism as transformations on the canonical variables generated through the 
Poisson brackets by 

Q[F]= J (fxUF{a,P,s) . (3.21) 
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This phase space functional is obtained from the charge (2.24) by using definition 
(3.2) for n. Note that the combination jiUi from equation (3.7) is invariant under 
these transformations, so from equation (3.9) n is invariant as well. It follows 
that Q[F] has vanishing Poisson brackets with the Hamiltonian and momentum 
constraints (3.11), confirming that Q[F] indeed generates a symmetry of the theory. 
Also observe that the symmetry generators (3.21) close under the Poisson brackets 
according to 

{Q[Fi],Q[F2]} = Q[F] , (3.22) 
where F = {dFi/da^){dF2/d(3^) - {dFi/d(3^){dF2/da^). 

3.3 Initial and boundary value problems 

A perfect fluid with equation of state p(n, s) coupled to the gravitational field 
is described by the canonical action (3.10) plus the canonical action for gravity. The 
Cauchy data for this system consist of the fiuid variables H, 9, (lis), a^, {11(3 a), 
and the canonical gravitational variables. These initial data cannot be specified 
independently, but must satisfy the Hamiltonian and momentum constraints. Any 
set of initial data that does satisfy these constraints can be transformed into an 
equivalent set by the symmetries (2.18), which are generated by the phase space 
functional Q[F] of equation (3.21). 

According to the analysis of section 2.2, the fields (p and 9 can be brought 
to zero on any spacelike hypersurface by a symmetry transformation. Thus, there 
is no loss of generality in choosing (p and 9 to be zero on the initial hypersurface. 
Moreover, the Lagrangian coordinates can be chosen to coincide with the coor- 
dinates a;* on the initial surface, so that = 6f. With these choices, equation 
(3.7) shows that the spatial components Ui of the covariant fluid four-velocity on 
the initial hypersurface are 

f^Ui = -Pi . (3.23) 

This reveals the geometrical signiflcance of the flelds Pa- with the choices (p = 9 = 0, 
a^^i = Sf allowed by symmetry, the fluid space covector components Pa are just 
—fi times the spatial components of the fluid four-velocity. 

Recall the definition Va := fiUa for the Taub current vector from equation 
(2.15). The result (3.23) shows that in specifying initial data, —Pa can be identi- 
fied with the spatial components Vi of the Taub vector. Thus, a complete set of 
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initial data for a perfect fluid consists of the (Eulerian) particle number density, 
the (Eulerian) entropy density, and the spatial part of the Taub vector, along with 
(p = 6 = and = (spatial coordinates). These initial data are then evolved 
according to the Hamiltonian differential equations of motion (3.19). 

Now consider the boundary value problem. Assume the spacetime manifold 
admits closed spacelike hypersurfaces so the boundary data is specified only on ini- 
tial and final hypersurfaces. One possible set of boundary data consists in specifying 
the canonical coordinates (p, 6, and on the initial and final hypersurfaces. These 
boundary data include 10 x oo^ boundary values, 5 x oo^ on the initial surface and 
5 X oo^ on the final surface, where oo^ is the number of space points. With these 
boundary data, the 10 x oo^ Hamiltonian first order differential equations of motion 
(3.19) generically determine a solution for the ten canonical field variables. 

Another possible set of boundary data consists in specifying the canonical 
momentum 11 along with the coordinates $ and on the initial and final hyper- 
surfaces. In this case, the data on the initial and final surfaces are related through 
the conserved charges Q[F] of equation (3.21). In particular, if and 11 are speci- 
fied initially and is specified finally, then the conserved charges Q[F] can be used 
to compute 11 on the final surface by considering a complete set of functions F{a"^). 
This means that the independent boundary data consist of only 9 x oo^ boundary 
values. With these boundary data, the 10 x oo^ Hamiltonian first order differential 
equations of motion generically determine the ten canonical field variables to within 
the symmetry transformation (2.19). In particular, the field (f is obtained only to 
within an additive constant along each of the oo^ flow lines. Other sets of boundary 
data are restricted by the conserved charges Q[F] as well. 

The distinctions among the various types of boundary data can be clarified by 
a simple example, namely, the free nonrelativistic particle. If the initial and final 
positions of the particle are given as boundary data, the equations of motion can 
be solved uniquely for the particle position as a function of time. But the initial 
and final momenta cannot be specified independently, because space translation 
invariance implies that the momentum is conserved. By specifying equal values for 
the intial and final momenta, the equations of motion can be solved to within a 
constant spatial translation of the particle. 

4. Velocity— potential representation 
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In equation (2.8a), namely 



Ua = -{<P,a + S9,a + Pa^^ ,a) / fJ' , 



(4.1) 



the fluid four-velocity is written in terms of various scalar flelds and their gradients. 
Expressions of this type are common in the literature on fluid dynamics. They are 
often called velocity-potential representations or Clebsch [29] representations of Ua, 
and the scalar flelds themselves arc called velocity potentials or Clebsch potentials. 
Two related questions naturally arise: Is the velocity-potential representation (4.1) 
sufficiently general to allow for any four-velocity, or does it restrict the four-velocity 
in some way? If the representation (4.1) is sufficiently general, then is it overly 
general in the sense that fewer potentials would be adequate? 

The first of these questions can be answered in the affirmative by explicitly 
constructing a set of velocity potentials 9, a"^, and (3a that correspond to an 
arbitrary timelike, unit normalized vector field t/", along with arbitrary spacetimc 
scalar fields s, T, and fi 0. Begin by defining the thermasy through the rela- 
tionship T = O^aU"^. More precisely, assign arbitrary values for 9 on some spacelike 
hypersurface, then define the value of 9 at any other spacetime point x by inte- 
grating the scalar T along a fiow line of f/" from the hypersurface to the point x. 
Likewise, choose arbitrary values for ip on some spacelike hypersurface then extend 
ip off this hypersurface by integrating the relationship f := (i — Ts = (p,aU"' along 
the flow lines of U"'. Now observe that the vector 



is orthogonal to the flow lines of U"'; that is, PaU"" = 0. Thus, Pa can be expressed as 
a linear combination of three independent basis vector fields, where the basis vectors 
span the subspace of spacetime vector flelds that are othogonal to U"'. Such basis 
vectors are the gradients of three scalar flelds a^, where are a set of Lagrangian 
coordinates that (uniquely) label the flow lines of t/". Therefore Pa can be written 
as Pa = PAo:^,a where Pa are scalar flelds. Inserting this expression for Pa into 
equation (4.2) yields the velocity-potential representation (4.1). 

The construction described above shows that the fluid four-velocity Ua always 
can be expressed in the velocity-potential representation (4.1) for any entropy per 
particle s, any nonzero chemical potential /x, and any temperature T such that the 



Pa-.^- 



{jJiUa + ip,a + S9^a) 



(4.2) 
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thermasy satisfies T = 9^aU°'. This construction also shows the extent to which the 
velocity potentials are arbitrary. In particular, (f and 6 can be chosen arbitrarily on 
a single spacelike hypersurface, and the fields can be chosen as any set of unique 
fiow line labels. These ambiguities in the velocity-potential representation appear 
as the invariances (2.19-21) of the perfect fluid action. 

Regarding the second question, it is indeed possible to reduce the number of 
potentials that appear in the velocity-potential representation by replacing the three 
pairs of fields a^, /?^, by a single pair (see references [4, 5, 15]). The demonstration 
invokes Pfaff's theorem [37] and the observation from section 2.2 that Pa can be 
viewed as a covector or one-form on the three-dimensional fluid space. According 
to Pfaff's theorem, can be written in terms of three fluid space scalar fields as 
Pa = <P,A + PoL^A- Then the velocity-potential representation (4.1) becomes 

Ua = -{<f,a + se^a + /^«,a)//i , (4.3) 

where the definition (f> := ^-\-(p has been used. Comparing equations (4.1) and (4.3) 
shows that the seven potentials a^. Pa, have been replaced by three potentials 
(fi, a, p. In effect, this reduction amounts to restricting the fluid space indices on 
and Pa to a single value. 

The perfect fluid action (2.5) admits the symmetry transformations (2.18) re- 
gardless of the range of values assumed by the fluid space indices. However, the 
identiflcation of as Lagrangian coordinates was used repeatedly in the interpreta- 
tion of those symmetries and the corresponding conserved charges (2.24). It should 
be emphasized that the reduction in the number of velocity potentials precludes 
this possibility of identifying as Lagrangian coordinates for the fluid. Thus, 
with just one pair of variables in place of and Pa, the physical interpretation of 
the symmetries and conserved charges is lost. Also recall that with three pairs a^, 
Pa, the flelds Pa can be given a direct geometrical interpretation by relating them 
to the spatial components of —Va to within a symmetry transformation (see section 
3.3). This allows for the speciflcation of initial data in terms of simple physical 
quantities, namely the Eulerian number density, the Eulerian entropy density, and 
the spatial part of the Taub current. 

Another disadvantage of the reduced representation (4.3) is that the velocity 
potentials (f, a, P are not always single valued functions on spacelike hypersurfaces. 
(This has been recognized in the context of nonrelativistic perfect fluids in reference 
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[38].) For example, consider a fluid whose four-velocity is described by the repre- 
sentation (4.1) with Pi = —o?^ (32 = a^, = 0. Such a fluid has nonzero vorticity 
with axis of rotation (in the isentropic case) in the direction e°'''''^(p^bCi^ ,cOi'^ ,d- Using 
the Pfaff reduction, P^ can be written as P^ = <fi,A + Pol^a where (p = R'^cj), a = R'^, 
and P = —(f). Here, R and (f) are polar coordinates in the a^-a^ plane of the fluid 
space, defined by R^ = (q;^)^ + (q!^)^ and tan0 = o? ja^ . Since the value of jumps 
by 27r along a spacelike loop surrounding the vorticity axis, the velocity potentials 
(p (and hence <^) and P are not single valued functions on space. 

This example shows that the representation (4.3), with potentials a, P that 
are single valued, does not allow for an arbitrary four-velocity. As a particular 
consequence, it can be shown [38, 39] that the fluid helicity is restricted to vanish. 
(Helicity is deflned, for example, in references [34] and [40].) On the other hand, the 
velocity-potential representation (4.1) can be used to represent any four-velocity 
with potentials that are single valued on all spacelike hypersurfaces. Since several 
coordinate charts may be needed to cover the fluid space with Lagrangian coordi- 
nates CK^, a more precise statement is that any four-velocity can be represented as 
in equation (4.1) by velocity potentials that are single valued within open subsets of 
spacelike hypersurfaces, where the open subsets contain flow lines from a single 
fluid space coordinate chart. 

In order to verify this claim, recall the discussion at the beginning of this 
section in which an arbitrary four-velocity Ua was constructed from the velocity- 
potential representation (4.1). The goal now is to argue that the potentials used in 
that construction are single valued. Assuming and d are chosen to be single valued 
on one spacelike hypersurface, then integration of the equations / = ^p^aU"' and 
T = d^aU"" yields spacetime flelds and 9 that are single valued on any spacelike 
hypersurface. Of course, s is always single valued since it has a direct physical 
interpretation. The Lagrangian coordinates from a given coordinate chart on 
the fluid space are also single valued functions on any spacelike hypersurface. From 
equation (4.2), the fields Pa on a hypersurface with coordinates are now defined 

by 

Pa = -Oi\ illUi + (p,i + Sd,i) , (4.4) 

where a\ is the inverse of a^^^. It follows that Pa arc single valued functions on 
any hypersurface since each field appearing on the right-hand-side of equation (4.4), 
including the spatial components Ui of the fluid four- velocity, are single valued. The 
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conclusion is that within a subset of any spacehke hypersurface that contains flow 
lines from a given fluid space coordinate chart, the velocity-potentials appearing in 
the representation (4.1) can be assumed to be single valued without any restriction 
on Ua. 

Another consequence of reducing the number of velocity potentials is that the 
counting of degrees of freedom (per space point) is changed. With all three pairs a^, 
(3a, the perfect fluid is described by flve pairs of unconstrained canonical variables 
and therefore has the expected flve degrees of freedom. Three of these correspond 
to the fluid's freedom of motion in three-dimensional space, while the remaining 
two degrees of freedom correspond to the fluid number density and entropy density. 
With just one pair of potentials replacing and (3 a, the system apparently has 
just three degrees of freedom. Evidently, this difference in the number of degrees of 
freedom arises because in the steps leading to the reduced representation (6.3), the 
equations of motion = and ^aU"" = were used in order to interpret 

Pa as a fluid space covector. The signiflcance of using the equations of motion can 
be understood by an analogy: Consider a nonrelativistic particle in three spatial 
dimensions moving in a central potential. The particle has three degrees of freedom, 
corresponding to motion in each of the three spatial directions. But the equations 
of motion show that the particle actually moves in a two-dimensional plane. Using 
this result, the problem can be reduced to that of a particle moving in its orbital 
plane, and the number of degrees of freedom is reduced to two. 

Finally, observe that if the term (p^a had been omitted from the representation 
(4.1), the Pfaff reduction would nevertheless lead to the representation (4.3) with 
an apparently trivial change of notation (p ^ 0. This reasoning has incorrectly 
lead to the conclusion (see the appendix of reference [5]) that the potential (p is 
unnecessary, and that a term ip,aJ°' is not needed for a valid action principle if 
all three pairs of flelds ck'^, /3a are used. The conclusion is not valid because (p 
is a function on the fluid space, and satisfles (p,aU°' = 0. This contradicts the 
result obtained by contracting equation (4.3) (with (p ^ p) with C7", which shows 
that ip,aU^ equals the chemical free energy /. Another (invalid) justiflcation [5] 
for dropping (p is that the variation of ip in the action just yields J",o = 0, and 
particle number conservation is already implied by the constancy of the labels 
along the fluid flow lines. The error in logic appears to lie in the interpretation 
of CK^ as particle labels, in which case their existence would imply conservation of 
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particle number. However, true particle labels are not continuous spacetime fields, 
so CK^ do not label individual particles. Rather, label the fluid flow lines that 
are physically determined by the average particle motion. The particle number 
conservation equation J'^^a = must be imposed as a separate equation of motion 
to insure that the number of particles within a flow tube (defined by a bundle of 
flow lines) is conserved. 

5. Action S with equation of state p(n, s) 
5. 1 Action and equations of motion 

Express the densitized fluid number flux vector as 

= -v^6"^'='^r7i23(«) a\ba\ca\d , (5.1) 

where 77123 is a function of a. The signiflcance of 77123(0:) can be seen by using the 
number flux vector nW^ = J"/ at a spacetime point x to construct a differential 
three-form on the spacelike hypersurface orthogonal to the fluid flow line at x. This 
three-form is 

^nWeabcddx^ A dx" A dx'^ = r]i2zda^ A da^ A da^ , (5.2) 

and is interpreted as the number of particles in the infinitesimal three-volume 
{l/3l)U"'eabcddx^ A dx'^ A dx'^. Then equation (5.2) shows that ?7i23(o) is the com- 
ponent of a three-form r] = {l/3l)r]ABcda'^ A da^ A da^ on the fluid space whose 
integral over a region B gives the number of fiuid particles whose fiow lines are 
included in B : 

rj — (number of particles in B) . (5-3) 

Under changes of coordinates (2.21a) in the fluid space, 77123 transforms as the 
component of a three-form. Correspondingly, expression (5.1) for J" is independent 
of the choice of fluid space coordinates. Note that locally, 77123 can be set to unity 
by an appropriate choice of coordinates. 

The equations of motion (2.8b,f), expressing conservation of particle number 
and constancy of the Lagrangian coordinates along the fluid flow lines, are im- 
mediately satisfled by virtue of the ansatz (5.1) for J". If the entropy per particle s 
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is given as a function of a^, the entropy exchange constraint (2.8c) is also automat- 
ically satisfied. Observe that sij is the fluid-space three-form whose integral over a 
region B gives the total entropy contained in the flow lines included in B. 

Given a three-form 77, a function s, and coordinates on the fluid space, 
a perfect fluid action S can be constructed that is a functional of the Lagrangian 
coordinates only. First consider the action S' = S — J d'^x{(pJ°' + sdJ°')^a, 
that differs from S of equation (2.5) by boundary terms. This action S' yields 
the equations of motion (2.8) that correctly describe a relativistic perfect fluid 
with equation of state p(n, s). From the action S", construct the functional S 
by substituting expression (5.1) for J"- and the function s{a) for s. The terms 
—'^J"',a — d{sJ"'),a + /3ACt^,aJ"' drop out, leaving the result 



where s = s{a) and J" is defined as a function of through equation (5.1). The 
procedure used to obtain this functional (5.4) does not guarantee its validity as an 
action principle. In particular, the equations of motion that follow by varying S 
with respect to a"^ are not necessarily equivalent to the original equations of motion 
from S or S' with J" expressed as a function of a'* through equation (5.1). Such 
a discrepancy would occur if the variations in induced by variations in are 
not the most general possible variations consistent with the constraint equations 
J"^a = 0, (sJ"),a = 0, and a"* qJ" = 0. However, in this case the representation 
(5.1) for J" is general, and the functional (5.4) does constitute a valid perfect 
fiuid action. This can be confirmed explicitly by varying the functional (5.4) and 
showing that the correct equations of motion follow. (Alternatively, the validity of 
the functional (5.4) as an action would be assured if it could be demonstrated that 
expression (5.1) and s — s{a) arise from solving the equations of motion (2.8a-d,f) 
for J", </?, 9, s and (3^ in terms of a^.) 

In varying the functional (5.4), observe first that variations with respect to the 
metric gab are unchanged by the substitution (5.1) for J". So S yields the correct 
stress-energy tensor (2.1). Variation with respect to gives 




(5.4) 




(5.5) 
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where 



abed 



o SAC D , I ^Vbcd \ b c d s: 

Sri^cnda ,bQ; ,0" ,d+ [ ]a ,hOL ,d^OL 



abed 



'i{riACDSoL^),bOi'^ ,cOi'^ ,d 



+ I 3 ^ „ I a ,bOi ,ea ,dda 



da' 



(5.6) 



The last term above is proportional to the four-form ^'r][Bcu,A] — 'nBCD,A — 'iVAicD,B], 
and must vanish because the fluid space is three-dimensional. Inserting 5J'^{a) into 
the variation (5.5) and integrating by parts yields the equation of motion 
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e''''"'Va.,bVAcna'',ca'',d-nTs,^ , 



(5.7) 



where Va = fiUa is the Taub current. Equation (5.2) implies rjAcoOt^ ,a(^'^ ,c(^° ,d 
nU'^ebacd, so the equation of motion simplifies to 



= ^e'^f'^^Ve-jU^aed - Ts,a 



(5.8) 



This is the combination (2.14) of the Euler equation and the entropy conservation 
equation. This analysis shows that the functional S is indeed a valid perfect fluid 
action: S has the correct stress tensor, incorporates particle number conservation 
and entropy conservation by virtue of expressions (5.1) for J" and s = s{a), and 
yields the Euler equation from its equations of motion. 

5.2 Hamiltonian form of the action S 
The momenta conjugate to are 



1 



(5.9a) 
(5.96) 



Here, is the matrix inverse of a^^i and 11 := J* is the spatially densitized Eulerian 
particle number density (3.9): 



n = \/lin{l + Uih'Wjf 



/2 



(5.10) 
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In principle, the next step in deriving the Hamiltonian is to solve equations (5.9) 
for as functions of Pa and . This is not possible for arbitrary equations 
of state p{n,s). The situation here is closely analogous to the expression of the 
Hamiltonian for the action S, discussed in section 3. In that case the Hamiltonian 
(3.11-12) is given as a function of the Lagrangian number density n and the spatial 
components Ui of the fluid velocity. In turn, n and Ui are determined as functions 
of the canonical variables through equations (3.7) and (3.9). For the case at hand, 
the Hamiltonian is again expressed by equations (3.11) and (3.12), namely 

= J d^x\^N^/h{p{l + WUi)+pWUi) + N'{-iinUi)^ . (5.11) 

This result is dictated by the fact that the fluid contributions to the Hamiltonian and 
momentum constraints are just the energy and momentum densities of the fluid and 
are given by appropriate projections of the fluid stress-energy-momentum tensor 
[17]. The quantities n and Ui that appear in the Hamiltonian (5.11) are determined 
as functions of the canonical variables a^, through equations (5.9b) and (5.10). 
For this purpose observe that, according to expression (5.1), 

n = ^^e^^SABc«^,i«"j«'',fc . (5.12) 

Thus, H only depends on a"^ and their spatial derivatives, not on a"^. Also note 
that the fluid contribution to the momentum constraint equals PaOl"^ ,i and is the 
canonical generator of spatial diffeomorphisms for the scalar fields and their 
conjugates Pa- 

The canonical equations of motion follow from the variation (3.18) of the fluid 
contribution to the momentum constraint, where H is the function (5.12) of a"^. 
This calculation shows that once again is given by equation (3.19c), and that 
Pa is given by 

- NVhTns,A + Ua\ ( (^^p^^^) ^ • (5-13) 

Equation (5.13) is just the space-time decomposition of the Lagrangian equation of 
motion (5.7). Also note that the flrst two terms in Pa above coincide with (11(3 a)' 
from equation (3.19f). 
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As mentioned in the introduction, the canonical Hamiltonian formulation of 
perfect fluids derived here is related to the Lie-Poisson Hamiltonian formulation 
by a Lagrangian to Eulerian map (see [18] and references therein). The result 
of this mapping is a Hamiltonian description of perfect fluids in which H, s, and 
the fluid momentum density are the fundamental variables. Since the Lie-Poisson 
brackets are not canonical, there is no corresponding action functional of the form 



For fixed fluid space tensors ry and s, the local coordinate expressions s{a) 
and r] 123 (ct) da^ A da'^ A da^ can be constructed from any set of coordinates 
on the fluid space. Under a change of coordinates Sa"^ = ^^{a), the functions 
^123(0;) and s{a) at a given value of (at a given coordinate location) change 
according to —(£^77) 123(0;) and —{£^s){a), respectively, and correspondingly the 
tensors 77123(0) da^ Ada'^ Ada^ and s{a) at a given point on the fluid space manifold 
remain unchanged. (Here, £^ is the Lie derivative in the fluid space along the vector 
fleld ^.) In this sense the functional S, which depends on a"^ only through the 
combinations 77123(0) da^ A da^ A da^ and s(o), is invariant under changes of fluid 
space coordinates. Note however that this invariance of S involves a transformation 
of the functions 77123(0) and s(o) at a given value of o'^, as well as a transformation of 
the field variables o"*. Consequently this invariance does not in general correspond 
to a conserved Nocthcr current. On the other hand, the subset of fiuid space 
coordinate transformations Sa"^ = that satisfy £^r] = and £^s = do not 
involve a transformation of the functions 77123(0) and s(o) at a given value of o"* 
and do give rise to conserved Noether currents. This same conclusion can be reached 
by considering the changes in 77123(0) da^ Ada'^ Ada^ and s(o) induced by a general 
field transformation da^ = where 77123(0) and s(o) are treated as fixed functions 
of o"^. Those changes are given by 



and show that 5" is invariant under transformations 6a^ = that satisfy £^r] = 
and £^s — 0. The Noether current associated with such transformations is obtained 



5.3 Symmetries of the action S 




{£^'n)i23{c()da^ Ada^ Ada^ , (5.14a) 
{£cs){a) , (5.146) 
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from a general variation of the action, which is 



SS = (terms giving the equations of motion) 



(5.15) 



For the transformations 5a^ = that leave S invariant, the last integral in ex- 
pression (5.15) must vanish when the equations of motion hold and the conserved 
current is 



The Noethcr charge is the integral of the time component of this current over a 
spacelike hypcrsurface. Using the definition (5.9) for the conjugate momenta, this 
charge can be written as 



Q[$] is independent of the spacelike hypcrsurface S for any fluid space vector ^ that 
leaves the tensors rj and s invariant under Lie transport. 

6. Other action functionals 

6.1 Equation of state p{fi, s) 

A relativistic perfect fluid with equation of state p(n, s) is described by his- 
tories that extremize the action (2.5), along with the gravitational action, under 
independent variations of J", (p, 9, s, a"^, and Qab- This choice of variables is 
convenient, but not unique. For example, replace J"'/\/^-g by nU^ with n and C/" 
varied separately, subject to the normalization condition Lf^Ua = — 1. This changes 
the individual equations of motion, but as a set the new equations are equivalent 
to the original equations (2.8). In particular, the stress tensor has the standard 
perfect fluid form (2.1) only when certain other equations of motion hold. 

The action (2.5) with J"'/\/—g replaced by nil"' is 



where n is an independent variable. Now, according to the flrst law (1.3), the 
chemical potential // is the function of n and s deflned by = dp/dn. If this 




(5.16) 




(5.17) 




(6.1) 
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equation could be inverted for n as a function of /i and s, then /i could be used as 
an independent variable in the action (6.1) by substituting n{fi, s) in place of n. In 
effect, this inversion is accomplished by using the pressure p as a function of fi and 
s to generate the functions n{ii,s), p{fi,s) through the first law in the form (1.4). 
That is, let p{fj,, s) specify the equation of state and determine n and p through 

dp 

n := - , (6.2a) 

dp 

p:=H—-p. (6.26) 

Inserting these relationships into the action (6.1) yields an action S(^p) for a perfect 
fluid with equation of state p{n, s). 

For the action 5'(p), it is convenient to use the Taub vector V"" = ^U"" as the 
independent variable in place of p, and f/". With this deflnition, the action with 
equation of state p{p, s) reads 

^(P) = J d^xV^!^p - (1^) {\V\ - y"(^,„ + se,a + f3.a\a)/\V\) I , (6.3) 



where p = p{p, s) and fj. \V\ = y—V^-gabV''- The equations of motion obtained 
by varying 9, s, a"*, and (3a yield equations (2.8b-f) respectively. The equation 
of motion obtained by varying V"' is 



The projection of this equation orthogonal to the fluid flow lines shows that {(p^a + 
sd,a + PAOi^,a) is proportional to V^, and then projection along the flow lines gives 

Va = -{^,a + se,a+pAa^,a) ■ (6.5) 

This is the equation of motion (2.8a), so the action S'(p) indeed yields the complete 
set of fluid matter equations (2.8). As shown in section 2, these equations of motion 
imply particle number conservation, conservation of entropy along the flow lines, 
and the Euler equation. The stress tensor derived from S'(p) has the perfect fluid 
form (2.1) when equation (6.5) holds. 
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The equation of motion (6.5) can be used to eliminate V"- from the action S'(p) 
yielding an action functional that is essentially the one found by Schutz [15]: 

^(P) = J d'^xV^Pil^, s) . (6.6) 

Here, /i is treated as the function of 9, s, a"*, (3a, and Qab determined by /U^ = 
—V°'Va, where Va is given by equation (6.5). The only difference between the action 
(6.6) and Schutz's action is a difference in the number of pairs of fields a"^ and 
Schutz uses just one pair, in accordance with the Pfaff reduction discussed in section 
4. The canonical form of Schutz's action is constructed in references [19] and [20]. 

6.2 Equation of state a{n, T) 

An action S(a) for a perfect fluid with equation of state specified by a(n, T), the 
physical free energy as a function of number density and temperature, is obtained 
by eliminating the entropy per particle s from the action (2.5). Thus, consider 
solving the equation of motion (2.8d) for s as a function of n and T = 9^aU°', then 
eliminating s from the action (2.5) by substituting the result s{n,T). The term 
—\f^p + sO^aJ"' becomes a function of n = \ J\/y/^ and T = 9^aJ°'/\J\ that is 
identified with no, where a is the physical free energy (1.2b). The resulting action 
for a perfect fluid with equation of state a(n, T) is 

Sia, = - j d^x[\J\a{\J\/^ , e,aJV\J\) - J%^,a + ^AOt^ ,a)) ■ (6.7) 

S^a) is a functional of J", 6*, ct'*, /S^, and gab- 

It is straightforward to show that the equation of motion 5S(^a)/^J°' = is 
equivalent to equation (2.8a), and the equations of motion obtained by varying 
9, a'*, and (3^ are just equations (2.8b,c,e,f). Of course, the equation of motion 
(2.8d) is missing, since S^a) does not depend on s. But with the identification of 
thermasy as in (2.9), this equation simply reiterates the relationship dictated by 
the first law of thermodynamics and is, in this sense, superfiuous. Therefore the 
equations of motion derived from the action S^a), along with the interpretation of 
the variables in a manner that is consistent with the first law (1.5), are complete 
in the sense that they imply particle number conservation, conservation of entropy 
along the fiow lines, and the Euler equation. In addition, the stress tensor obtained 
from the action (6.7) has the perfect fiuid form (2.1). 
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The action discussed by Kijowski et al. [14] can be obtained from the action 
S(^a) — J d'^x{(pJ°')^a by substituting expression (5.1) for J'*. The terms —(pJ°'^a + 
I^aOl^ ,aJ°' drop out, leaving 

S,a, = - j d^x[\J\ a(| J|/v^, e,aJV\J\)] ■ (6.8) 

Variations with respect to the metric gab and the thermasy 9 are unchanged by this 
substitution, so S^^a) yields the correct stress tensor (2.1) and the conservation of 
entropy (2.8c). A calculation similar to the one appropriate for S shows that vari- 
ations of S'(a) with respect to the Lagrangian coordinates yield equation (2.14), 
which implies the Euler equation. The Hamiltonian form of S^a) and the symmetries 
and conserved charges can be found along the same lines as the analysis found in 
section 5 for S. The Hamiltonian form of this action also has been considered in 
reference [41]. 

6.3 Isentropic fluids and dust 

Isentropic fluids are perfect fluids with a constant entropy per particle s. The 
first law of thermodynamics in the form (1.5) indicates that isentropic fluids are 
described by an equation of state of the form 

a(n,T) = ^-sT, (6.9) 
n 

where s is the constant value of the entropy per particle. Inserting this equation 
of state into the action (6.7) yields a functional like S of equation (2.5), but with 
two differences: the function p only depends on n = \ J\/y/—g, and s appears as a 
fixed constant, not a variable. Thus, s9 can be absorbed into <f by the change of 
variables + s6, which reduces the combination </? ^ + s9^a that appears in 

the action to (p'^^- This leads to the isetropic fiuid action 

-^isentropic = j d^x[-^ p{\J\/ ^) + if[,J'' + P^a"" ^aJ''} , (6.10) 

which is a functional of J", p' , a"^, /3a, and gab- In terms of the new variable p>' , the 
equation of motion (2.8a) obtained by varying becomes = pUa + p>' a + /^ACi'^.a- 
Contracting with t/" gives p — 'p[a^°' ^^'^ shows that ip' is a "potential" for the 
chemical potential p. 
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Dust is a particular case of an isentropic fluid in which the energy density 
p is proportional to the number density n and the pressure (2.7) is zero. The 
proportionality constant is the rest mass-energy of a fluid particle, which equals 
the chemical potential jj. = p/n. Inserting the equation of state p = p,n into the 
isentropic fluid action (6.10) yields 

Sdust = J d^x[-f,\J\ + <p,aJ'' + pAa^aJ"-} , (6.11) 

where the prime has been dropped from (p. An alternative to this dust action is ob- 
tained by replacing \ J\/y/—g in the action (6.11) with the function [{\J\/\/--g)'^/n + 
n] /2 and treating n as a new dynamical variable. This is justified because the action 
(6.11) is recovered when n is eliminated from the new action by using the solution 
n = \J\/^/^ to the equation of motion for n. After introducing the variable n into 
the dust action, J'* can be eliminated by using the J°' equation of motion. This 
yields an action for dust that is a functional of n, (p, a^, and (3 a, namely 

^dust = / d^xV^niUag'^'Ut + 1) , (6.12) 

where Ua := -(v,a + (^a^^ ,a) / IJ- 

The Hamiltonian form of the dust action (6.11) or (6.12) is just the canonical 
action (3.10) but with the terms containing 9 and s omitted: 

^dust = j dtd^x^ip + (n/3^)a^ - NW'^'' - N'nT^'"^ ■ (6.13) 

The dust contributions to the Hamiltonian and momentum constraints can be found 
explicitly, because the chemical potential appearing in equation (3.7) for Ui has no 
n dependence. Thus, from equations (3.14a) and (3.15), 

Hr*=nv?,. + (n/?^)a-^,, , (6.14a) 

W'^'' = [{iiuf + nt^'h'^nf^'] ^''^ . (6.146) 

In reference [26], the action (6.12) and its canonical form (6.13) are used to analyze 
the quantum theory of gravity coupled to dust. 

The equation of state (6.9) also can be inserted into the action functional 
(6.8). This yields an isentropic fluid action 

-^isentropic = - / d^ x\s/^ p(\J\l %f^) \ , (6.15) 
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where a total divergence (boundary term) has been discarded. Here, J"" is a function 
of the Lagrangian coordinates as given by equation (5.1). Speciahzed to a 
pressureless dust equation of state, this action becomes ^dust = — J d^x iJ,\J\. In 
the canonical form of this dust action, the dust contribution to the momentum 
constraint is P^a^^j, and the dust contribution to the Hamiltonian constraint is 
given by equation (6.14b) where 11 is a function of through equation (5.12). 
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Appendix 

The symmetries of the action (2.5) are transformations among the potentials 
ip, s, 9, Pa, and that leave the one-form —jjU = dip + sdd + (S^da^ and the 
entropy per particle s invariant at each spacetime point. That is, the action satisfies 
S[ipjSj6j(3A,Oi^] = S[ip' J s' ,9' , P'^, a"^'] if the primed and unprimed variables are 
related by 

dip + sd9 + P^da^ = dip' + s'd9' + p'^da^' , (A.la) 
s = s' . {A.lh) 

In general, the primed variables can be nonlocal functions of the unprimed vari- 
ables. Thus, for example, ip'{x) can depend on the unprimed variables at spacetime 
points other than x. The analysis that follows is restricted to the case of ultralocal 
transformations in which the primed variables at x depend only on the unprimed 
variables at x. In particular, the primed variables are not allowed to depend on the 
derivatives of the unprimed variables. 

Equation (A. la) can be viewed as the expression of a canonical transformation 
with generating function ip — ip' for a fictitious phase space with coordinates 9, a"^ 
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and momenta s, /J^. In terms of the true phase space variables discussed in section 
3, the symmetry transformations (A.l) are given by 



Udif + {Us)de + {Ul3^)da^ = U'dif' + {UsYdO' + {U/3^yda 

n = n' , 
(n.) = (n.)' . 



A' 



{A.2a) 
{A.2b) 
{A.2c) 



Equation (A. 2a) expresses a Mathieu transformation among the canonical fluid vari- 
ables at each space point. (Mathieu transformations are canonical transformations 
with zero generating function; that is, they preserve the form Pidq'^ [1].) The extra 
restrictions (A. 2b) and (A. 2c) can be analyzed by computing Poisson brackets of 
the primed variables with one another using the definition of the Poisson brackets 
in terms of the unprimed variables. For example, condition (A. 2b) along with the 
Poisson bracket relationship {(p',Il'} = 1 lead to d(p'/d(p = 1; therefore. 



where / is independent of (p. The Poisson bracket of (p' and (lis)' along with 
condition (A. 2c) show that / is independent of 9. A similar analysis gives 



where g, like /, is a function only of 11, (lis), a^, and {11/3 a)- It also follows from 
conditions (A. 2b) and (A. 2c) that a' and (11/3^)' are independent of (p and 9. 
Using the results (A. 3) and (A. 4), equation (A. 2a) can be written as 



where is defined by IIjF = 11/ + {Ils)g + (n/3^)'Q;^'. This relationship can be 
investigated by treating various combinations of primed and unprimed variables as 
independent. For example, assume that (11/3^)' as a function of 11, (lis), a"^, and 
(11/9^) can be inverted for (11/3^) as a function of 11, (lis), ct^, and (11/3^)'. Then 11, 
(lis), q;'*, and {11/3 a)' can be chosen as independent variables and equation (A. 5) 
immediately yields 



^' = ^ + f , 



{A.3) 



9' = 9 + g, 



{AA) 



d(njP) = fdU + gd{Us) + (n/3^)da^ + a^'d(n/3^)' , 



(A.5) 




{A.6a) 
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(HA) = , (A.6c) 

These equations imply 

= u{dj^/du) + {Us) {dj^/d{Us)) + (up^) {dj^/d{up^y) , {a.7) 

and show that is homogeneous of degree zero in the momenta H, (lis), and 
{UP^y. As a consequence, the functional dependence of can be expressed as 
= :F(q!'*, (ns)/n, (n/3^)'/n) and equation (A. 6a) can be written as 

Other choices of independent variables lead to alternative sets of transformations, 
in the same way as canonical transformations are categorized into various types. 
The infinitesimal symmetry transformations are obtained by setting 

= {a^ {U(3^y) /U + eF , (A.9) 

where e is an infinitesimal parameter. Equations (A. 8) and (A.6b,c,d) then become 

dF OF 

where F has functional dependence F{a'^, (lis)/!!, (n/?^)/!!). The factors of 11 can 
be eliminated from these relationships, leading directly to the expression (2.18) of 
the infinitesimal symmetry transformations in terms of the fiuid potentials </?, s, 9, 
(3a, and a^. 
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